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An analysis of irreversibilities generated due to combustion in an adiabatic combustor burning wood was 
conducted. This was done for a reactant mixture varying from a rich to a lean mixture. A non-adiabatic 
non-premixed combustion model of a numerical code was used to simulate the combustion process 
where the solid fuel was modelled by using the ultimate analysis data. The entropy generation rates due 
to the combustion and frictional pressure drop processes were computed to eventually arrive at the 
irreversibilities generated. It was found that the entropy generation rate due to frictional pressure drop 
was negligible when compared to that due to combustion. It was also found that a minimum in irre¬ 
versibilities generated was achieved when the Air-Fuel mass ratio was 4.9, which corresponds to an 
equivalence ratio of 1.64, which are lower than the respective Air-Fuel mass ratio and equivalence ratio 
for complete combustion with theoretical amount of air of 8.02 and 1. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Incomplete combustion process of solid fuels such as those 
found in circulating fluidised bed combustors is difficult to model 
as compared to complete combustion [1 ]. The reason for this is that 
for complete combustion the proportion of each combustion 
product can easily be predicted analytically through the stoichio¬ 
metric equation [1—4], whereas this cannot be the case for 
incomplete combustion. This makes it hard to predict the perfor¬ 
mance of a combustion unit as measured against desired output. 
The AF mass ratio (Air-Fuel) can be selected such that combustion 
occurs, but for incomplete combustion the proportions of the 
combustion products cannot be predicted analytically as 
mentioned earlier. For an incomplete combustion process the 
combustion products mixture will contain unburnt fuel in the form 
of char and maybe volatile gases released from the devolatilisation 
of the solid fuel. There is negligible molecular oxygen in the com¬ 
bustion products mixture [1—4], In the case of complete combus¬ 
tion with excess air there will be molecular oxygen in the 
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combustion mixture, but there will not be any fuel components 
that resulted from the devolatilisation process [1 ], assuming proper 
mixing occurred. Because the combustor is adiabatic the total 
enthalpy across it does not change, but entropy increases [ 1 ]. A way 
to alleviate this predicament is to conduct experiments or run 
numerical simulations that predict the proportion of each com¬ 
bustion product. Elfasakhany et al. [5] undertook experimental and 
numerical studies to find reasons for emissions of unburned vola¬ 
tiles and particles of pulverised wood in combustion products 
mixture, and they found that the emissions were sensitive to fuel 
particle size and combustor configuration. Different AF ratios will 
results in different exergies or availability of the energy released as 
measured by means of a combustion products average tempera¬ 
ture. Chen et al. [6] compared the gasification of raw biomass, 
torrefied biomass and coal using pure oxygen, and found that there 
is an optimum oxygen-fuel ratio for each fuel at which gasification 
is most efficient. Tran et al. [7] found that torrefaction of stump 
biomass resulted in increased heating values and grind-ability of 
the fuel when compared to untreated biomass, but torrefaction 
made the treated fuel less reactive. Zhang et al. [8] experimented on 
gasification of caking coal in a pilot fluidised bed gasifier, where 
they found that the optimal proportion of oxygen used for pre¬ 
oxidation of 37% of the total oxygen supply greatly enhanced the 
gasification rate. The availability of the energy released from the 
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Nomenclature 

Sgen 

entropy generation rate [W/K] 



S h 

viscous dissipation energy source term 



S m 

solid phase to gas phase mass transfer source term 

Abbreviation 

T 

temperature [K] 

CFD 

computational fluid dynamics 

u 

speed [m/s] 



~v 

velocity [m/s] 

Symbol 


X 

molar fraction 

a 

velocity constant term [s -1 ] 

Y 

mass fraction 

A 

cross-sectional area [mm 2 ] 



AF 

Air-Fuel mass ratio 

Greek symbol 

c 

molar specific heat [kj/kmol K] 

P 

density [kg/m 3 ] 

c 

specific heat [kj/kg K] 

e 

void fraction 

Cp 

specific heat at constant pressure [kj/kg K] 

e 

turbulent dissipation rate 

C 

atomic carbon 

it 

one of the mixture fraction variance transport equation 

Cd 

one of the mixture fraction variance transport equation 


constant 


constant 

P 

molecular viscosity [Pa s] 

Cg 

one of the mixture fraction variance transport equation 

Pt 

turbulent viscosity [Pa s] 


constant 

T 

stress tensor 

CO 

carbon monoxide 



C0 2 

carbon dioxide 

Subscript 

/ 

mixture fraction 

0 

superficial 

f 12 

mixture fraction variance 

Air 

air 

F 

body force 

char 

char 

g 

gravitational acceleration magnitude [m/s 2 ] 

CO 

carbon monoxide 

~g 

gravitational acceleration [m/s 2 ] 

C0 2 

carbon dioxide 

h 

specific enthalpy [kj/kmol] 

f 

fuel 

h° 

enthalpy of formation [kj/kmol] 

g 

gas 

H 

atomic hydrogen 

gen,min 

minimum entropy generation rate 

H 

height [mm] 

gen(fi) 

entropy generation rate due to heat 

H 

total enthalpy [kj] 

gen(p) 

entropy generation rate due to frictional pressure drop 

h 2 

molecular hydrogen 

h 2 

molecular hydrogen 

h 2 o 

water vapour 

h 2 o 

water vapour 

T 

unit vector 

i 

counting index 

/ 

irreversibility [W] 

l 

maximum number of products 

k 

turbulent generate rate 

m 

maximum number of reactants 

L 

length [mm] 

m 

mixture 

m 

mass flow rate [kg/s] 

mf 

minimum fluidisation 

M 

molecular weight [kg/kmol] 

n 2 

Molecular nitrogen 

n 

stoichiometric coefficient 

0 2 

molecular oxygen 

N 

atomic nitrogen 

P 

product species 

N 

entropy generation number 

r 

reactant species 

n 2 

molecular nitrogen 

ref 

environment reference 

0 

atomic oxygen 

s 

solid species 

0 2 

molecular oxygen 

sd 

dense zone solid species 

p 

pressure [Pa] 

se 

exit zone solid species 

R 

universal gas constant [kj/kmol K] 

t 

total of solid and gas 

S 

specific entropy [kj/kmol K] 

th 

theoretical amount of air ratio 

s° 

absolute entropy [kj/kmol I<] 




combustion process can be quantified through the computation of 
an entropy generation rate as given by Bejan [9], Le Roux, Bello- 
Ochende and Meyer [10] used the technique of total entropy gen¬ 
eration minimisation to find the optimum operating point of a 
turbine, and to optimise the geometric parameters of a recuperator 
and receiver of an open and direct solar thermal Brayton cycle. 
Mwesigye, Bello-Ochende and Meyer [11] numerically investigated 
the entropy generation rate due to heat transfer and fluid friction in 
a parabolic trough receiver. They found that a minimum entropy 
generate rate existed for every combination of concentration ratio 
and inlet temperature, which corresponds to an optimum flow rate. 
Le Roux, Bello-Ochende and Meyer [12] used the dynamic trajec¬ 
tory optimisation method to find optimised geometric parameters, 
minimum irreversibility rates and maximum receiver surface 


temperatures for an open and direct solar thermal Brayton cycle. 
For an adiabatic circulating fluidised bed combustor the entropy 
generation rate is due to frictional pressure loss across the 
combustor and the combustion process. The frictional pressure 
drop across a conduit for a fluidised system as given by Kunii and 
Levenspiel [13] remains the same once minimum fluidisation has 
been attained. Gupta and Reddy [14] analysed the effect of pressure 
on the bed-to-wall heat transfer in a pressurised circulating flui¬ 
dised bed riser by modifying a cluster renewal model developed for 
atmospheric circulating fluidised bed riser columns. The model 
they developed agreed well with experimental results and litera¬ 
ture. Han et al. [15] found that pressurisation improved the capture 
of CO2 when using CaO in a fluidised bed gasifier of sawdust. 
Aisyah, Ashman and Kwong [16] found that coal fly-ash oxygen 
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carriers were effective in converting carbon dioxide into carbon 
monoxide when used in chemical looping combustion of synthesis 
gas. Aghbalou, Badia and Ilia [17] conducted an exergetic optimi¬ 
sation of a solar collector and thermal storage system where they 
used a phase change material. The optimisation resulted in an 
analytical model that compared well to experimental data they 
generated a priori. Kousksou et al. [18] conducted a second law 
analysis of latent thermal storage for a solar system where they 
studied the behaviour of phase change material in the form of its 
melting temperature. Exergy analysis of phase change materials 
was also conducted by Li et al. [19] where they developed a 
mathematical model for the thermal storage of energy from a solar 
collector system. They developed the model based on the finite¬ 
time thermodynamics method. A review of exergy based perfor¬ 
mance evaluation of phase change materials latent heat thermal 
storage systems was performed by Jegadheeswaran, Pohekar and 
Kousksou [20] where they looked at the influence of design pa¬ 
rameters and operating conditions on optimisation of the thermal 
storage systems. Brandvoll and Bolland [21 ] used exergy analysis on 
results from simulations to show that a chemical looping com¬ 
bustion process has lower irreversibilities than a conventional 
combustion process with air. Srinivas et al. [22] performed energy 
and exergy analyses of pressurised and atmospheric circulating 
fluidised bed char combustion units where they identified the 
pressure ratio and maximum value of the gas turbine inlet tem¬ 
perature as the parameters that have a dominant impact on the 
combined cycle plant performance. The exergy efficiency of a blast 
furnace process was found to have high exergy efficiency, but the 
compressing and preheating of the blast was found to be charac¬ 
terised by high exergy losses due to the associated combustion 
process in a study conducted by Ziebik and Stanek [23], Ziebik and 
Stanek [24] also looked at the thermal improvement of a blast 
furnace process through the increase of blast temperature and its 
oxygen enrichment. They did this through energy and exergy an¬ 
alyses. Nur Izyan and Shuhaimi [25] used exergy analysis on a crude 
oil distillation unit to find strategies to reduce fossil fuel depletion 
problem, they found that the first inlet furnace had the highest 
exergy loss and that the best strategies are to reduce the heat loss 
from furnace stack and an overall cleaning schedule of crude pre¬ 
heat train. Gupta and Reddy [14] performed energy and exergy 
analyses on a combined-cycle power plant all in the aim of mod¬ 
ernising the plant. They found that the combustion chamber, gas 
turbines and heat recovery steam generators were the main sources 
of irreversibilities. The irreversibilities of Proton Exchange Mem¬ 
brane fuel cells were analysed by Cihan, Hacihafizoglu and Kahveci 
[26] as a function of operating pressure, pressures of the electrodes, 
current density, etc. Bapat, Majali and Ravindranath [27] did an 
exergetic evaluation of two quintuple effect evaporators with and 
without heat recovery, and found that an increase in exhaust inlet 
temperature results in deterioration in exergetic performance of 
the evaporator with heat recovery and the quality of end product. 
Ghazikhani et al. [28] used a double pipe heat exchanger with 
counter current flow to recover exergy by cooling the exhaust from 
a diesel engine. They found that the recovered exergy increased 
with increasing load and engine speed. They also found that when 
they used the recovered exergy, the brake specific fuel consump¬ 
tion decreased markedly. Ghazikhani et al. [29] investigated the 
improvement in performance and reduction in emissions due to 
ethanol additives in a two stroke SI (spark ignition) engine. They 
found that there is a reduction in brake specific fuel consumption 
and there is a reduction in pollutants emitted, with a 35% reduction 
in CO emission being the highest. In a related paper Ghazikhani, 
Hatami and Safari [30] found that for 5% added ethanol into the fuel 
the second law efficiency was higher for 2500 and 3000 revolutions 
per minute settings. Chintala and Subramanian [31 ] used an exergy 


analysis on a hydrogen-diesel engine to determine the maximum 
available work and irreversibility due to mixing, combustion, un¬ 
burned and friction. They found that the maximum available work 
increased with H2 addition due to reduction in irreversibility due to 
combustion. They also found that the addition of H2 had no effect 
on the irreversibility due to friction and mixing, whereas the irre¬ 
versibility due to unburned fuel decreased due to the high tem¬ 
perature combustion that result with the addition of H 2 . Jafarmadar 
[32] found that even though the exergy efficiency decreased when 
they increased the mass fraction of exhaust gas recirculation in an 
indirect injection diesel engine from 0% to 30%, the cumulative 
irreversibility in the combustion chamber decreased. Ting, Hung 
and Guo [33] investigated the entropy generated by nanofluid with 
streamwise conduction in microchannels. They found that there 
existed minimal entropy generation at certain range of the low- 
Peclet number flow regime. The need to combat rise in CO and 
CO2 emissions from coal combustion has put pressure on utilities 
and governments of countries where the main source of electricity 
is from coal combustion. Substitutes for coal combustion that can 
result in lower emissions of CO and CO2 are some of the mitigation 
solution under consideration at the moment. These substitutes can 
be complete replacement of coal with wood as a source of heat in a 
boiler or a blending of wood with coal. However combustors usu¬ 
ally do not completely burn the solid fuel, especially circulating 
fluidised combustors. In order to investigate the viability of wood as 
a source of heat in a typical combustor setting, the analysis of ir¬ 
reversibilities generated during the combustion of pitch pine in an 
adiabatic circulating fluidised bed combustor was conducted as the 
Air—Fuel mass ratio was varied to achieve a reactant mixture 
varying from a rich mixture to a lean mixture. This was achieved by 
fixing the fuel mass flow rate at one value, but we reduced or 
increased the air mass flow rate. 

2. Mathematical model 

The study was of an adiabatic combustor burning pitch pine 
wood as illustrated in Fig. 1. The adiabatic combustor consisted of a 
cylinder with an adiabatic wall, inlets for the solid fuel and primary 
air inlet for the air used to devolatilise the solid fuel and secondary 
air inlets if needed to complete the combustion of the fuel. The 
different states of complete/incomplete combustion were attained 
by changing the mass flow rate of the incoming oxidising air, 
thereby varying the AF. The chemical formula of the solid fuel as 
computed in Ref. [1 ] is given by C4.9H7.2O2. 

The lower heating value of the fuel from Ref. [ 1 ] was given as 
16,091.2 kj/kg and the molecular weight of the fuel was calculated 
as 98 kg/kmol. The specific heat at constant pressure was also given 
as 1.680 kj/kg K. 

The availability of the energy released from the combustion 
process was analysed for scenarios on both sides of the theoretical 
air complete combustion point, i.e. the incomplete combustion 
with air deficiency and complete combustion with excess air. 
Complete combustion of pitch pine wood occurs at an AF of 8.02 as 
was determined in Ref. [1[. The availability was analysed as a 
function of AF. 

2.1. Exergy analysis 

The irreversibility generation rate that results from the com¬ 
bustion and frictional pressure drop processes taking place in an 
adiabatic combustor is as given by Eq. (1) 

^-|(Enr i (h-r r efl) ri -Enp i (S-T ref s) pj ) (1) 
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Flue gas 



Fig. 1. Schematic of an adiabatic combustor (adopted from Ref. [1]). 


Since the enthalpy terms cancel out in Eq. (1 ), the irreversibility 
generation rate can also be defined, by Eq. (2) 

I = TWSgen (2) 


2.1.1. Entropy generation rate 

Entropy generation rate is a term that can be used to quantify 
the inefficiencies of thermodynamic systems as they convert en¬ 
ergy from one form to another. Complete elimination of this term is 
not possible for real systems; however its minimisation could result 
in maximisation of the exergy available from these systems and 
hence improved efficiencies. The entropy generation rate term due 
to combustion was defined [9] in terms of the change in entropy of 
the system due to reactions. On the other hand the entropy gen¬ 
eration rate due to frictional pressure drop was defined in terms of 
the pressure drop of the system [9], Each of the two terms was 
defined and derived in turn as presented below. 


2.1X1. Friction pressure drop entropy generation rate. The friction 
pressure drop entropy generation rate term is made up of the 
contribution due to the air-fuel mixture as it flows through the 
riser. The contribution due to the air-fuel mixture flow is given by 
Eq. (6). 


Sgen(p) 


m t AP 

PmTm 


(6) 


AP is the pressure drop across the combustor, rh t is the total air and 
fuel mass flow rate, p m is the mixture density, T m is the mixture 
temperature. The pressure drop across the riser due to the air-fuel 
mixture [13] is given by Eq. (7). 

AP = Lmf (l - e mf ) (p f - P g )g (7) 


e m f is the minimum fluidisation solid volume fraction, pf is the char 
density, p g is the gaseous density, L m f is the minimum fluidisation 
height and g is acceleration due to gravity. The pressure drop can 
also be expressed in terms of the total height of the riser by using an 
expression found in Ref. [13] for an air-solid mixture undergoing 
fast fluidisation inside a vertical column as given by Eq. (8). 


where S gen is the entropy generation rate. The specific entropy is / \ £ sd — e se / ,\ 

given by Eq. (3). ^mf ~ «mf ) = - - -+ H te sd - H { (e sd -e s ) (8) 


Spi = 5pi(P)-Pln(P*p,/P re f) (3) 

The static pressure in the combustor is the same as the reference 
static pressure, i.e. P = P re f, and the molar fraction of each product 
species, x P j, is given by Eq. (4) 



n P j is the stoichiometric coefficient of combustion product species i. 
The same expressions apply for reactants, h is the specific enthalpy, 
Spj is the absolute entropy, R is the universal gas constant, p is the 
static pressure, P re f is the reference static pressure of the environ¬ 
ment and T re f is the reference temperature of the environment. The 
specific entropy of the solid fuel was computed by defining [1 ] it as 
given by Eq. (5). 

S = cln(7>/298) (5) 

c is the molar specific heat of the solid fuel. 

The entropy generation rate S gen is due to the combustion and 
the pressure drop phenomena. Generally the entropy generation 
rate term is not defined in detail in terms of its constituents. In 
Section 2.1.1 we give the detailed expressions of the entropy gen¬ 
eration rate term constituents. 


e s d is the solid volume fraction in the bed region of the riser, £ se is 
the solid volume fraction at the exit of the riser, a is an empirical 
velocity term, H t is the total height of the riser, Hf is the freeboard 
height and £* is the minimum value that £ Se can assume. Kunii and 
Levenspiel [13] also give an expression for Hf as given by Eq. (9). 

H f w. [ln(e sd - 4) - ln(e se - 4)] /a (9) 


When Eqs. (7)—(9) are substituted into Eq. (6), the pressure drop 
is eliminated and results in an expression as given by Eq. (10). 


_ nt t (pf-pg)g[&L%e + H t £ sd + ln(|^ 4 ) (^V^)] 
Sgen(p) - ~ 7 . \... ' 

[CsePf + (1 — e S e)PgJ T m 


( 10 ) 


The mixture density was computed as expressed by Eq. (11). 
Pm = esePf + (1 ~ e se )p g (11) 

The empirical velocity term a can be expressed in terms of the 
fuel mass flow rate and the AF. Firstly the empirical velocity term is 
given by a = 10/u 0 , where u 0 is the superficial velocity. But the 
superficial velocity is given asuo = (1 - £ S d)u g , where u g is the gas 
velocity. The gas velocity in turn is expressed as Eq. (12) 


u g = thgjpgA 


(12) 
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Given the fact that the total mass flow rate is expressed by Eq. 
(13), and that the Air—Fuel mass ratio is expressed as AF = m s /rhf, 
the superficial velocity can be written as expressed by Eq. (14). 

nit = m g + rhf (13) 


U 0 = (1 - e sd ) 


m f AF 

"JgA 


(14) 


Hence the empirical velocity term can be expressed in terms of 
the fuel mass flow rate and AF as given by Eq. (15). 


2.1.2. Entropy generation number and relative entropy generation 
rate 

The entropy generation rate at any equivalence ratio was ana¬ 
lysed further by means of two ratios. An equivalence ratio of 1 
means complete combustion with theoretical amount of air, 
whereas an equivalence ratio above 1 signifies a rich mixture and 
that below 1 signifies a lean mixture. The first ratio is the entropy 
generation number, N s [9] which is the quotient of the entropy 
generation rate at any equivalence ratio and the minimum entropy 
generation rate, and is as given by Eq. (21). 


10p g A 

(1 - e sd )m f AF 


(15) 



(21) 


When Eq. (15) is substituted into Eq. (10), the entropy genera¬ 
tion rate due frictional pressure can be expressed in terms of the 
fuel mass flow rate and AF as given by Eq. (16). 


Since S gen min is the best possible solution at which the most 
exergy can be generated, N s measures the penalty paid in terms of 
lost exergy when deviating far from the point where minimum 


Sgen (p) — 


m f (AF +1) ( Pf - P g )g[ ^- e i ( o fi - e / )miAF +H t e sa + ln(^|) (^)(1 - esd )m f AF] 

(csePf + (1 — e se)PgJ T m 


(16) 


2.1.12. Combustion entropy generation rate. The entropy genera¬ 
tion rate due to the combustion process is given by Eq. (17). 


_ rhfAs 
sen(h > ~ M f 


(17) 


The change in entropy S gen ^ accounts for the increase in en¬ 
tropy as the fuel is oxidised by the oxygen in the air. The change in 
entropy is constituted by contributions from gaseous species and 
solid species, i.e. As = As g + As s . The contribution to the entropy 
change from gaseous species is given by Eq. (18). 


As g = E n Pl( s ° - Rln (*))pj - rj(s° -Mn(x)) rJ 


(18) 


entropy generate rate occurs. The second ratio is relative generation 
rate, which is the quotient of the entropy generation rate at any 
equivalence ratio and the entropy generation rate when combus¬ 
tion occurred with theoretical air ratio, and is as expressed by Eq. 

( 22 ). 

Nr = (22) 

S genth measures the entropy generation rate when complete 
combustion occurs with theoretical amount of air. This point also 
happens to be where maximum possible temperature can be ach¬ 
ieved. So N r can be a measure of gained or lost efficiency by 
sacrificing attainment of the maximum possible temperature in 
order to maximise exergy generation. 


Mj is the molecular weight of species j. The contribution to the 
entropy change from solid species is as expressed by Eq. (19). 

ASs = ^ n charCcharM char In (V char / 298 ) 

-|n fCf M f ln^gg) (19) 

c c har is the specific heat constant for the char particles and Cf is the 
specific heat constant for the solid fuel. When Eqs. (18) and (19) are 
substituted into Eq. (17) the entropy generation rate due to com¬ 
bustion can be expressed as given by Eq. (20). 


2.2. CFD (computational fluid dynamics) model 

The adiabatic combustor produced 241 kW t h when complete 
combustion of pitch pine wood occurs with theoretical amount of 
air. The combustor had a diameter of 300 mm and a height of 
7000 mm as illustrated in Fig. 2. 

The mass flow rate of the incoming solid fuel was fixed at a value 
of 0.015 kg/s so as to have a common maximum heat generation 
rate of 241 kW t h since the lower heating value is fixed and unique 
for any particular fuel. ANSYS Fluent 14 [34] was used to simulate 
the combustion process inside the adiabatic combustor. The com¬ 
bustion process was modelled using the non-premixed combustion 
model because the solid fuel is only specified in ultimate analysis 


mj 

M f 


' E %(s“ — Rln(x)) p . — E n rj (s° 
7=1 7 = 1 

^-| nf C f Mf ln (rf /298 ) 


^ln(X)) rJ - + ^ n char c char^char l n ^char/298^ 


(20) 











J. Baloyi et al. / Energy 70 (2014) 653-663 


Combustion Products Outlet 



Fig. 2. Schematic of the adiabatic combustor showing boundary placements. 

data as given in Table 1, and the combustion model gives an option 
of specifying the fuel as an empirical stream. The inlet diffusion 
option was selected. 

The continuity and momentum equations [34] as expressed by 
Eqs. (23) and (25) respectively were solved for under steady state 


condition. 

V-(pT?7) = v-(gvj)+s m (23) 

V- (W 2 ) = V- (gvp) + C g fi t (v/) 2 - C d p JF (24) 

V-(/>m?) SvteVp + V-(t) + pt + F (25) 

The stress tensor [29] is given by Eq. (27). 

f « I *[(W + W T ) - Jv-T??] (26) 


p is the molecular viscosity of the continuous phase, F is the 
interactive body forces between the dispersed and the continuous 
phases, S m is the source term accounting for the mass transfer from 
the solid phase to the gas phase and i is a unit vector, f and f ' 2 are 
the mixture fraction and its variance, and are computed by applying 


O' Pressure outlet 


Adiabatic wall 


Secondary air mass inlets 




Fig. 3. Schematic of the geometry used in ANSYS Fluent 14 to model the adiabatic 
combustor. 

an assumed shape probability density function (/3-function) when 
modelling the turbulence-chemistry interaction. 

The energy for the non-premixed combustion model [34] was as 
expressed by Eq. (27). 

V.(plTH) = yQW)+S h (27) 

The total enthalpy H is given by Eq. (28). 

m T 

H=E Y J f ( C Pj dT + h !( T refj)) (28) 

J ' = 1 Tr¬ 


iable 1 

Pitch pine ultimate analysis data (adapted from Ref. [35]). 
Element Value 


Table 2 

Mesh independence analysis results. 


H 

0 


LHV 

Molecular weight 


59.0 (%) 

1.2 (%) 

32.7 (%) 

1.13 (%) 

24,220 (kj/kg) 
16,091.2 (kj/kg) 
98 (kg/kmol) 


Number of cells S ge n (W/K) A (%) 


438,192 362.6218 

609,857 362.6756 0.01 

655,783 362.7186 0.01 

747,055 364.1874 0.40 

1,188,339 361.6429 -0.70 
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Yj is the mass fraction, c p j is the specific heat and constant pressure 
and hj>(r ref j) is the enthalpy of formation of the jth species. Sh is the 
source term due to viscous dissipation. The use of the chosen 
combustion model required the use of a turbulent model because 
the combustion model is a mixture model. To this end the k—e 
turbulent model with enhanced wall function option was chosen 
for all simulations. 

The adiabatic combustor had four mass flow inlets for the solid 
fuel as shown in Fig. 3, each with equal mass flow rate of 
0.00375 kg/s. The inlet temperature of the fuel at each inlet was set 
at 600 K. The temperature value was selected to be same as the 
devolatilisation temperature of pitch pine [35], The oxidising air 
was split into the primary and secondary air, each making up half of 
the total air used in the oxidation of the fuel. The reason for split¬ 
ting the oxidising air was to enable the modelling of incomplete 
combustion to take place. The primary air entered the combustor at 
the base with a single fixed mass inlet, with half of the total air mass 
flow rate. The secondary air four fixed mass inlets were situated a 
quarter of the height up the combustor, each having an eighth of 
the total air mass flow rate. The inlet temperature at all eight air 
inlets was set at 298 K. The wall of the combustor was adiabatic. 
The data that was extracted from the simulations were the tem¬ 
perature at the pressure outlet boundary and the combustion 
products molar fractions. Fig. 3 shows the placement and types of 
boundaries applied on the combustor. 

Mesh independence of the solution was carried out for different 
mesh cell numbers. The Figure of merit used to quantify mesh 



Fig. 5. The entropy generation rate profile as a function of AF in an adiabatic 
combustor. 


independence was A, which was used to quantify the change in 
entropy generation rate as the number of cells in the computation 
domain was increased. The entropy generation rate was computed 
for the case of complete combustion with theoretical amount of air. 
It was computed as A — (S gen ,• — S geru _i )/S g en,i-i, and the results 
are tabulated in Table 2. 

The mesh used for the numerical simulations had 747,055 un¬ 
structured polyhedral cells. The Presto scheme was used to solve for 
pressure and second-order upwind schemes were used to solve for 
the continuity, energy, turbulence and mass fraction. The results 
obtained when using Eqs. (16) and (20) to process data from CFD 
simulations and data from thermodynamics tables data [2] were 
compared for the case of complete combustion with varying 
amount of excess air, that is lean mixtures. It was found that the 
results from CFD data compared well with results from thermo¬ 
dynamics tables data as shown in Fig. 4. 

3. Results and discussion 

The entropy generation rate was almost entirely due to the 
combustion process, with negligible contribution from the 
component from frictional pressure drop as can be seen from Fig. 5. 

The change in the entropy generation rate as a function of AF 
was plotted as shown in Fig. 6. The figure shows that the entropy 
generation rate reaches a minimum at an AF of 4.9. This is well 
illustrated by the close up plot as shown by Fig. 7. This ratio occurs 



Fig. 7. A close up of the entropy generation rate profile as a function of AF in an 
adiabatic combustor. 
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Table 3 

Computed values of the entropy generation rate and irreversibilities. 

AF Equivalence S gen(p) (W/K) S gen{h) (W/K) S gen (W/K) / (W) 


4 2.005 

4.1 1.956098 

4.2 1.909524 

4.3 1.865116 

4.4 1.822727 

4.5 1.782222 

4.6 1.743478 

4.7 1.706383 

4.8 1.670833 

4.9 1.636735 

5 1.604 

5.1 1.572549 

5.2 1.542308 

5.3 1.513208 

5.4 1.485185 

5.5 1.458182 

6 1.336667 

6.5 1.233846 

7 1.145714 

7.5 1.069333 

8.02 1 

8.42 0.952494 

8.82 0.909297 

9.63 0.832814 

10.43 0.768936 

11.23 0.714159 

12.03 0.666667 


0.237390033 3644.42 
0.240265992 2824.117687 
0.243068457 2036.776795 
0.245796486 1544.100364 
0.248450379 1242.88416 


0.251030896 

0.253537704 

0.25597112 

0.258330312 

0.260615659 

0.262827439 

0.264965702 

0.26703045 

0.269021683 

0.270939401 

0.272783605 

0.280901916 

0.287182364 

0.291624947 

0.294229666 

0.294988965 

0.294220374 

0.292275546 

0.284734689 

0.27255259 

0.255665558 

0.234073594 


916.5697472 

653.7978266 

369.1511782 

274.588215 

265.0651819 

268.7854397 

272.828133 

277.0432976 

280.6182435 

284.5333752 

288.6333938 

304.1216374 

318.76 

332.9877197 

350.3486205 

364.1874005 

371.7057442 

378.5330356 

390.9738209 

403.2551109 

414.5166439 

424.663812 


3644.66 

2824.358 

2037.02 

1544.346 

1243.133 

916.8208 

654.0514 

369.4071 

274.8465 

265.3258 

269.0483 

273.0931 

277.3103 

280.8873 

284.8043 

288.9062 

304.4025 

319.05 

333.2793 

350.6429 

364.4824 

372 

378.8253 

391.2586 

403.5277 

414.7723 

424.8979 


1086.11 

841.6587 

607.0319 

460.2152 

370.4535 

273.2126 

194.9073 

110.0833 

81.90427 

79.06709 

80.17638 

81.38174 

82.63848 

83.70441 

84.87169 

86.09404 

90.71196 

95.08 

99.31724 

104.4916 

108.6158 

110.856 

112.8899 

116.595 

120.2512 

123.6021 

126.6196 


when there is incomplete combustion inside the adiabatic 
combustor as complete combustion with theoretical amount of air 
occurs at 8.02. Fig. 6 also shows the variation of the generated ir¬ 
reversibilities with AF. which were computed from the entropy 
generation rate as given by Eq. (2). 

The computed values of the entropy generation rate due to the 
combustion process, frictional press drop, the overall entropy 
generation rate and the irreversibilities generated were also tabu¬ 
lated in Table 3 for all the plotted AF and corresponding equivalence 
ratios. 

Fig. 8 also shows the variation of entropy generation rate, but as 
a function of equivalence ratio. An equivalence ratio of 1 means 
complete combustion with theoretical amount of air, whereas an 
equivalence ratio above 1 signifies a rich mixture. Fig. 8 also shows 
that the entropy generation rate reached a minimum at an equiv¬ 
alence ratio of 1.64 where incomplete combustion occurs. Both 
Figs. 6 and 8 show that when you have an increase in richness of the 





Fig. 8. Plot of entropy generation rate profile as a function of equivalence ratio in an 
adiabatic combustor. 



Fig. 9. Plot of entropy generation number and relative entropy generation rate as a 
function of AF in an adiabatic combustor. 


reactants mixture or deficiency in air the entropy generation rate 
drops off gradually until minimum is reached. Then the entropy 
generation rate rises rapidly as the AF is decreased gradually. 

The entropy generation rate was normalised with respect to the 
value computed when combustion occurred with theoretical 
amount of air Nr, and this plotted in Fig. 9. Fig. 9 shows that an AF of 
about 12 (which is a 50% excess air) resulted in a relative increase in 
entropy generation rate of 16.5%. However an AF of 4.5 resulted in 
about 2.5 fold increase in entropy generation rate. Fig. 9 also shows 
that at an AF of 4.7 entropy was generated at the same rate as when 
there is complete combustion with theoretical air ratio. A similar 
plot for the relative entropy generation rate is plotted in Fig. 10, but 
as a function of equivalence ratio. 

Figs. 9 and 10 also show that the minimum entropy generation 
rate is about 73% of the value when combustion occurred with 
theoretical air ratio. 

However the relative entropy generation rate was also 
computed with respect to the minimum value as an entropy gen¬ 
eration number, N s , and is plotted in Fig. 9. Fig. 9 shows that the 
entropy generated when there is complete combustion with 
theoretical air ratio of 8.02 is about 37% more than the minimum 
rate. Combustion at an AF of 4.5 resulted in an entropy generation 
rate that is almost 3.5 times that of the minimum rate which 
occurred at an AF of 4.9. The entropy generation number was also 
plotted as a function of the equivalence ratio as shown in Fig. 10. 

The temperature of the exiting combustion products mixture 
was plotted as a function of the AF as shown by Fig. 11. As expected 
the highest temperature is achieved when combustion occurred 



Fig. 10. Plot of entropy generation number and relative entropy generation rate as a 
function of equivalence ratio in an adiabatic combustor. 
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Fig. 11. Plot of combustion products mixture temperature as a function of AF in an 
adiabatic combustor. 


with theoretical air ratio of 8.02, and it decreases as either a lean or 
rich mixture is combusted. 

The variation of molar fraction of combustion species was plotted 
as a function of the AF as shown in Fig. 12. The molar fraction for N2 
was not included in the plot because it is far bigger than the others, 
hence it would obscure the variations in the other molar fractions, 
and the fact that nitrogen was assumed to be inert to the combustion 
process. Fig. 12 could hold a clue as why the minimum entropy 
generation rate occurred at an AF of 4.9. It can be seen from the 
figure that at about an AF of 4.9 the molar fractions of H 2 and H 2 0 are 
equal. This ratio also happened to be about the same point at which 
the molar fraction of C0 2 is at its lowest for a rich mixture. The other 
species are negligible in comparison. A consequent effect of this is 
that the combined contribution to total entropy generation rate 
from H 2 , H 2 0 and C0 2 is the lowest at an AF of 4.9 than at any other 
AF, since at an AF of 4.9 the molar fractions that are significant other 
than that of N 2 are those of C0 2 , H 2 0 and H 2 . This then result in the 
lowest entropy generation rate occurring at an AF of 4.9. As 
mentioned in Ref. [1] the combustion model implemented in the 
CFD software of [31] over-predicts the presence of 0 2 . It was also 
mentioned in Ref. [1] that this has the effect of lowering the com¬ 
bustion products mixture temperature due to the presence of more 
species available to share in carrying the heat energy. 

As mentioned earlier the entropy generation rate due to fric¬ 
tional pressure drop is negligible compared to that due to the 
combustion process. However for the geometric configuration used 
in this study the entropy generation rate due combustion cannot 
tell if the flow process is physically possible or not. For this we used 
the entropy generation rate due to frictional pressure drop to see 
where it became negative to know whether it would or not occur in 
the physical world. For the current geometric and flow 




Fig. 13. Plot of entropy generate rate due to frictional pressure drop as a function of AF 
in an adiabatic combustor. 

configuration, Fig. 13 shows that the flow process is physically 
possible because S gen (p) is positive throughout the AF range under 
consideration. Because it is negligible in magnitude in comparison 
to the entropy generation rate due to combustion, the entropy 
generation rate due to frictional pressure drop was plotted on its 
own as a function of the AF as shown in Fig. 13. Whereas a mini¬ 
mum is observed for the overall entropy generation rate, Fig. 13 
predicts a maximum for the frictional pressure drop entropy gen¬ 
eration rate. The maximum value occurred at an AF of about 8, 
which is almost an AF of 8.02, where complete combustion 
occurred with theoretical air ratio. Section plots of contours of 
temperature (in Kelvins) through the adiabatic combustor for the 
case with AF of 4.9 and 8.02 are plotted respectively in Figs. 14 and 
15. From the two figures it can be seen that the average adiabatic 
temperatures plotted in Fig. 11 were arrived at by averaging over a 
temperature range at the outlet of the combustor. 
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Fig. 12. Plot of variation 1 
adiabatic combustor. 


tion of combustion product species in; 


Fig. 14. Section plot of filled temperature contours (Kelvins) for the case with AF of 4.9. 
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4. Conclusions 

The entropy generation rates due to the combustion and fric¬ 
tional pressure drop processes were analysed in order to calculate 
the irreversibilities generated in the adiabatic combustor. It was 
found that the entropy generation rate due to frictional pressure 
drop was negligible when compared to that due to combustion. The 
analysis of the irreversibilities generated during the combustion 
process of an adiabatic combustor with a fixed geometry was 
conducted and it was found that a minimum was reached when the 
Air-Fuel mass ratio was 4.9. At about this Air-Fuel mass ratio the 
molecular hydrogen and water vapour were at equal amounts of 
molar fraction and the carbon dioxide molar fraction was at its 
minimum for a rich mixture. This resulted the combined contri¬ 
bution from H2, H2O and CO2 to the entropy generation rate being 
the lowest at an AF of 4.9 than at any other AF; hence a minimum 
entropy generation rate since they have the most significant molar 
fractions other than N2. 
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